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Abstract. Semiclassical periodic-orbit theory and closed-orbit theory represent a quantum spectrum as a 
superposition of contributions from individual classical orbits. Close to a bifurcation, these contributions 
diverge and have to be replaced with a uniform approximation. Its construction requires a normal form 
that provides a local description of the bifurcation scenario. Usually, the normal form is constructed in flat 
space. We present an example taken from the hydrogen atom in an electric field where the normal form 
must be chosen to be defined on a sphere instead of a Euclidean plane. In the example, the necessity to base 
the normal form on a topologically non-trivial configuration space reveals a subtle interplay between local 
and global aspects of the phase space structure. We show that a uniform approximation for a bifurcation 
scenario with non-trivial topology can be constructed using the established uniformization techniques. 
Semiclassical photo-absorption spectra of the hydrogen atom in an electric field are significantly improved 
when based on the extended uniform approximations. 

PACS. 03.65.Sq Semiclassical theories and applications - 32.60.+i Zeeman and Stark effects 



1 Introduction 

Periodic-orbit theory and, as a variant for atomic 

photo-absorption spectra, closed-orbit theory [30] have 
proved to be the most powerful tools available for a semi- 
classical interpretation of quantum mechanical spectra. 
Closed-orbit theory provides a semiclassical approxima- 
tion to the quantum response function 
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which involves the energy eigenvalues E n and the dipole 
matrix elements (i\D\n) connecting the excited states n to 
the initial state i. Semiclassically, the response function 
splits into a smooth part and an oscillatory part of the 
form 

g osc (E)=^Ac. . (E) e iS - W h , (2) 

CO. 

where the sum extends over all classical closed orbits start- 
ing from the nucleus and returning to it after having been 
deflected by the external fields. 5 C . . is the classical action 
of the closed orbit. The precise form of the recurrence am- 
plitude A c .o. depends on the geometry of the external field 
configuration, but it can in any case be calculated from 
classical properties of the closed orbit. Periodic-orbit the- 
ory, which can be applied to any quantum system possess- 
ing a classical counterpart, represents the quantum den- 



sity of states as a sum over the classical periodic orbits in 
complete analogy with (J2J. 

Both variants of the semiclassical theory share the 
problem that the semiclassical amplitudes A c .o. diverge 
when the corresponding classical orbits undergo a bifur- 
cation. This failure occurs because close to a bifurcation 
the classical orbits cannot be regarded as isolated in phase 
space, as the simple semiclassical expansion {2J assumes. 
It can be remedied if the individual contributions of the 
bifurcating orbits to the closed-orbit sum (J2J) are replaced 
with a uniform approximation that describes their collec- 
tive contribution [o]|d]|7||o| . For the generic codimension- 
one bifurcations of periodic orbits [^||7||^J as well as closed 
orbits |9tlU| . uniform approximations are available in the 
literature. Typically, bifurcations of higher codimensions 
are not encountered classically because they split into se- 
quences of codimension-one bifurcations upon an arbitrar- 
ily small perturbation. Semiclassically, however, these se- 
quences must be described collectively by a single uni- 
form approximation if the individual bifurcations are suf- 
ficiently close |llU12Ulrilll4| . It can thus be necessary in 
practice to construct uniform approximations for bifurca- 
tion scenarios that are much more complicated than the 
elementary bifurcations of codimension one. 

The crucial step in the derivation of a uniform approxi- 
mation is the choice of a normal form <l> a {x) de- 
scribing the bifurcation scenario. Its stationary points, i.e. 
the solutions of d<P a (x)/dx = 0, must be in one-to-one cor- 
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respondence with the classical closed orbits. In addition, 
as the parameter a is varied, the stationary points must 
collide and change from real to complex in accordance 
with the closed-orbit bifurcations. Once a suitable normal 
form is known, the uniform approximation assumes the 
form 

&{E)= 1(a) e iSo ^/ h , (3) 

where 

and the normal form parameter a, the action So and the 
unknown function p(x) are determined as functions of the 
energy E such that the uniform approximation asymptot- 
ically reproduces the closed-orbit formula J2J for isolated 
orbits if the distance from the bifurcation is large. 

Usually, the choices of a and So are unique. By con- 
trast, the asymptotic condition fixes the values of the am- 
plitude function p(x) only at the stationary points of <P a , 
so that there remains considerable freedom in its choice. 
Customarily, a low-order polynomial that interpolates be- 
tween the given values is chosen. This decision can be 
justified by observing that a uniform approximation is 
needed only in the vicinity of the bifurcation, where the 
stationary points of the normal form are close. Since, in 
the spirit of the stationary-phase approximation, only the 
neighbourhood of the stationary points of <P a contributes 
to the integral Q, a Taylor series approximation that is 
accurate in that region can be expected to be appropriate. 

So far, the configuration space of a uniform approxi- 
mation, i.e. the domain of the normal form variable x and 
the range of the integral (@J could always be chosen to be 
a one- or multidimensional Euclidean space. In this paper, 
we describe an example where this is impossible because 
for topological reasons the bifurcating orbits cannot be 
represented as stationary points of a function defined on 
a Euclidean space. This result is surprising because uni- 
form approximations, as well as the underlying normal 
forms, are needed to describe closely spaced orbits and bi- 
furcations, and in this sense are local constructions. The 
topology of the configuration space, on the other hand, re- 
flects its global characteristics. The necessity to base a uni- 
form approximation on a configuration space of a suitable 
topology thus involves a subtle interplay between local 
and global aspects of the classical phase space structure. 
We demonstrate that the established techniques for the 
construction of uniform approximations are, in principle, 
applicable to uniform approximations on a topologically 
non-trivial configuration space, but that additional com- 
plications arise when non-local aspects of the dynamics 
have to be dealt with. 

The bifurcation scenario to be discussed arises, e.g., 
in the hydrogen atom in an electric field. This is an inte- 
grable system, and photo-absorption spectra can in prin- 
ciple be obtained by the quantization of classical tori |151 
I16| . However, while the torus quantization is restricted to 
integrable systems, closed-orbit theory is intended to be 
applicable to arbitrary atomic systems exhibiting either 
regular, chaotic, or mixed classical behaviour. Therefore, 
it is of vital interest to see how it can be applied to this 



apparently simple example. It turns out that to a closed- 
orbit theory approach the hydrogen atom in an electric 
field poses the same difficulties in principle as a generic 
mixed system because the closed orbits undergo plenty of 
bifurcations. We thus start with a description of the closed 
orbits in that system in section|21 In section[3]we focus on 
the bifurcations at low energies and show that in this spec- 
tral region a uniform approximation is required that goes 
beyond the known uniform approximation for the Stark 
system 17,18,19,20.21 . This uniform approximation is 
constructed in section 0] where we also demonstrate that 
it must rely on a sphere instead of a Euclidean plane as 
its configuration space. It is compared with the previous 
uniform approximation, and in section[S]its merits are as- 
sessed by the calculation of semiclassical photo-absorption 
spectra. 

2 Closed orbits 

The bifurcations of closed orbits in the hydrogen atom in 
an electric field follow a simple scheme that was described 
by Gao and Delos [22- Complementary to that work, we 
demonstrated in a recent publication |21| that to a large 
extent the closed orbits and their parameters can be com- 
puted analytically. We will now summarize the essential 
features of the bifurcation scenario. For details, the reader 
is referred to [2Tj . 

If we assume the electric field F to be directed along 
the z-axis, the Hamiltonian describing the motion of the 
atomic electron reads, in atomic units, 

H=\ V 2 --+Fz, (5) 
2 r 

where r 2 = x 2 + y 2 + z 2 . It can be simplified by means of 
its scaling properties: If the coordinates and momenta are 
scaled according to x = F x / 2 x, p = F~ 1 / 4 p, the classical 
dynamics is found not to depend on the energy E and the 
electric field strength F separately, but only on the scaled 
energy E = F- 1 / 2 E. 

The qualitative simplicity of the bifurcation scenario 
is due to the well-known fact |23II24| that the Hamilto- 
nian Q is separable in semiparabolic coordinates p\, p2, 
ip, where 

pi = \/f - 1 , p 2 = \/f + I , (6) 

and ip is the azimuth angle around the electric field axis. 
The angular coordinate <p is ignorable, so that its conju- 
gate momentum, the angular momentum L z around the 
field axis, is conserved. For closed orbits, which pass through 
the nucleus, L z — 0. The dynamics of p\ and pi is then 
given by the equations of motion 

pf= p\ + 2Ep 2 + Cl , 
p' 2 = -pi + 2Ep 2 + c 2 , 

where E is the energy and the separation constants c\ and 
C2 describe the distribution of energy between the p\ and 
P2 degrees of freedom. They are restricted by 

ci + c 2 = 4 . (8) 
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For real orbits the separation constants are constrained to 
the interval 



< Cj < 4 . 



(9) 

If they fail to meet (0, either p\ or p 2 assumes complex 
values, so that a closed "ghost" orbit in the complexified 
phase space is obtained. 

In the simplest cases, one of the pj vanishes identically. 
The electron then moves along the electric field axis. If 
ci = 0, the downhill coordinate p± is zero. The electron 
leaves the nucleus in the uphill direction, i.e. in the di- 
rection of the electric field, until it is turned around by 
the joint action of the Coulomb and external fields and 
returns to the nucleus. After the first return, the uphill 
orbit is repeated periodically. As the external field pro- 
vides an arbitrarily high deflecting potential, the uphill 
orbit is closed at all energies. 

The second case of axial motion is obtained if ci = 
4, which corresponds to p 2 = and thus to a downhill 
motion opposite to the direction of the electric field. The 
downhill orbit closes only if the scaled energy E is less than 
the scaled saddle point energy E$ — —2. If E > Eg, the 
electron escapes to infinity and the orbit does not close. At 
scaled energies E < E§, the downhill orbit, as the uphill 
orbit, is repeated periodically. 

In general, the motion of the electron will not be di- 
rected along the electric field axis and will consist of a 
superposition of the downhill and uphill motions. In this 
case a closed orbit is obtained if the downhill and uphill 
motions close at the same time after fc oscillations in the 
downhill direction and I oscillations in the uphill direction. 
A closed orbit characterized by integer repetition number 
(fc, I) must satisfy [22 



fcp_ K-(mi) = lp + K(m2) 



with 



p+ = \/-E+ \JE 2 - Cl , mi = +4, 

P+ 



(10) 



(11) 



C-2 



P- = \j-E+yE 2 + c 2 , m 2 = --^ 

and fC{m) denoting the complete elliptic integral of the 
first kind [23 . 

For given repetition numbers fc, I and scaled energy 
E, if we use that C2 = 4 — ci by iJBJl, equation (fTTHl can 
be solved for the separation parameter ci that specifies 
the initial conditions of the orbit. A solution only exists 
if I > fc. Once ci is known, all orbital parameters rele- 
vant to closed-orbit theory can be found analytically. In 
particular, the classical action of the closed orbit is [21] 



S=-kpll(m 1 )--lp 3 _I(m 2 ) 



(12) 



with 



I(m) = (m - l)K,(m) + (m + l)£ (m) (13) 

and the complete elliptic integrals of the first and second 
kinds JC(m) and £(m) |25]. It is real for both real and 
ghost orbits. 



If the repetition numbers fc and / are fixed, for suffi- 
ciently low scaled energies E the value of the separation 
constant ci obtained from (|10|) is greater than 4, so that 
the (fc, ?)-orbit is a ghost. As E increases, ci decreases 
monotonically. At the energy E^ est where Ci = 4, the 
(fc, /^-orbit becomes real. It is generated (as a real orbit) 
in a bifurcation off the fc-th repetition of the downhill or- 
bit, which is characterized by ci = 4. As E is further 
increased, the orbit moves away from the downhill orbit 
and approaches the uphill orbit. At the energy -Edest where 
ci = 0, the orbit (fc, I) collides with the l-th repetition of 
the uphill orbit and is destroyed. At E > Edest the orbit 
is a ghost again. 

The bifurcation energies E gcn and Edest can be deter- 
mined from (|10fl if ci = 4 is prescribed for the generation 
of an orbit, or ci = for its destruction. For these cases, 
equation (|1C)|> simplifies to 



2 3/2 




(14) 



in terms of the dimcnsionless variable e = 4/E 2 — (E^/E) 2 
for the generation of the (fc, I) orbit and 



2 3/2 



JC 



(15) 



1 71^/1 + VTTl V(i + \/i~R / 

for its destruction. 

3 Bifurcations at low energies 

From the discussion of the closed orbit bifurcations in the 
preceding section it is evident that there is only a single el- 
ementary bifurcation scenario in the hydrogen atom in an 
electric field, viz. a non- axial orbit being born or destroyed 
in a collision with an axial orbit. Uniform semiclassical 
approximations describing this type of bifurcation have 
been derived by several authors |17U18U19II21| . However, 
they are only applicable if the individual bifurcations are 
sufficiently isolated. If they are close, more complicated 
uniform approximations that provide a collective descrip- 
tion of several bifurcations must be found. Thus, there are 
two spectral regions where the construction of uniform ap- 
proximations poses a particular challenge: 

First, as the Stark saddle energy E§ — —2 is ap- 
proached from below, the downhill orbit undergoes in- 
finitely many bifurcations in a finite energy interval before 
it is destroyed. The periods of the closed orbits created in 
this way, as the period of the downhill orbit itself, grow ar- 
bitrarily large. The actions, however, remain finite, so that 
the action differences between the orbits stemming from 
two successive bifurcations become arbitrarily small. Sim- 
ilar cascades of isochronous bifurcations have been found, 
e.g., in the diamagnetic Kepler problem |26| and in Henon- 
Heiles type potentials |27| . Their uniformization remains 
an open problem. A detailed discussion of this cascade in 
the Stark system must be referred to future work. 
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Fig. 1. Bifurcations of non-axial orbits off the downhill (+) 
and uphill (x) orbits in a scaled period vs. scaled energy plot. 
Bifurcation energies and scaled periods of the bifurcating orbits 
are indicated. Bifurcations of the downhill orbits only occur 
below the Stark saddle point energy at Es = — 2. Dashed lines 
indicate first, second and third multiples of T ■■ 



2E. 



Second, non-axial orbits with winding ratios z = l/k 
close to one are created and destroyed at low energies, 
and their generation and destruction are closely spaced 
(see figure 0. Therefore, at low energies a uniform ap- 
proximation describing both the creation of a non-axial 
orbit from the downhill orbit and its destruction at the 
uphill orbit is required. It will be derived in section As 
a prerequisite, to quantitatively assess its necessity, the 
asymptotic behaviour of closed-orbit bifurcations at low 
energies will be investigated in detail in this section. 

Bifurcation energies can be determined from equations 
<|14|) and (|15fl . As a winding ratio z = l/k close to one 
corresponds to low bifurcation energies and small values 
of the parameter e = (Es/E) 2 , approximate bifurcation 
energies can be obtained by expanding the right hand sides 
of l|14|l and (|15|) in a Taylor series around e = 0. We then 
find the generation and destruction energies 



E, 



V3 



35 



'gen 



Edc 



2^/z~T 16^3 

V3 23 
2^/z~l 16V3 



Vz=l + 0(z- 1) 3/2 
+ O (z - l) 3/2 



(16) 



for an orbit with a given rational z, so that the energy 
difference 



AE = E, 



dest 



29 
8%/3 



VI^l + O (z - if' 2 (17) 



indeed vanishes as z — > 1. 

However, the crucial parameter determining whether 
or not the bifurcations can be regarded as isolated is not 
the difference of the bifurcation energies, but rather the 
action differences between the non-axial orbit and the two 
axial orbits. For the simple uniform approximation to be 
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Fig. 2. Schematic illustration of the bifurcation scenario with 
the definition of the critical quantities. Thick solid line: non- 
axial real orbit, thick dashed line: non-axial ghost orbit, thin 
long-dashed line: downhill orbit, thin short-dashed line: uphill 
orbit. Vertical lines indicate the bifurcation energies. 



applicable, at any energy at least one of the two action 
differences must be large. In other words (see figure EJ, at 
the critical energy where the actions of the downhill and 
the uphill orbits are equal, the difference between their 
action and the action of the non-axial orbit must be large. 

To begin with, the actions of the downhill and uphill 
orbits are given by 



S 



3tt 



down 
up 



-5/2 



' >r>7f (-E)- 9 / 2 ±^^(-E)-^ 2 (18) 



+ 



32^2 
45045tt 

8192 V2 



512^ 
-E)~ 17 ' 2 +<D(-E)~ 19 > 2 . 



The critical energy for the orbit (fc, I) is characterized by 
the condition fcSdown — ^u P , so that S& own /S uv 
solve this equation for the critical energy 



z. We 



E r 



V3 V3 



2 X /F^T 
209^ 



69%/3 
64 



(z-1) 



3/2 



256 



{z-\fi 2 + o(z-iyi 2 



(19) 



As expected, the critical energy lies between the two bi- 
furcation energies. If it is substituted into H18[). the critical 
action of the axial orbits 



37rfc 



(z-1) 5 / 4 



4557rfc 
576 • 3V4 



4 • 3V4 
(z-l) 9 / 4 + 0(z-l) 13 / 4 



(20) 



is obtained. 
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To determine the critical action of the non- axial orbit, 
the separation parameter ci for this orbit at the critical 
energy must be calculated from equation l|10|) . which can 
be rewritten as 



p- /C(mi) 







p+/C(m 2 ) 

Substituting the critical energy (|19fl into 1|21[1 . we find 



(21) 



29(ci 



p_ /C(mi) _ 

p+ /C(m 2 ) ^ ~ 24 

219 c? - 876 ci + 526 



+ 



144 



(z-lf 



(z-lf 



17951 c? - 107986 c? + 188112 ci - 76688 4 



9216 



+ 0(z- l) 5 



(22) 



We now make a power series ansatz for ci , insert it into lll'l't 
and solve for the coefficients to find 

175 175 

ci = 2+ — {z-l)~ — (z-lf + 0(z-lf. (23) 

From the knowledge of the critical energy and the separa- 
tion parameter, the critical action of the non-axial orbit 
is found to be 



° non 3I/4 v* L > ^ A . ql/4 ^ ^ 



115l7Tfc 

576 • 3 x /4 



4 ■ 3V4 

(z - I) 9 / 4 + 0(z - 1) 



13/4 



(24) 



This result agrees with (|20() in leading and next-to-leading 
order. The critical action difference is 



AS C 



29 !l74 (* - 1) 5/4 + 0(z - I) 9 / 4 



24 • 3 1 / 4 



(25) 



with q — I — k. If the critical action difference is written 
as a function of the critical energy given by H19fl . it reads 



AScrit — 



29nq 
32\/2 



Merit)- 572 +0(-£crit)- 7/2 



(26) 



Note that due to delicate cancellations occurring during 
the calculation, all terms of the series expansions given 
above are needed to obtain the leading-order result (|2"6*|l . 
Due to these cancellations, the critical action difference 
is proportional to a high power of the energy, so that it 
vanishes fast at low energies. 

Figure shows numerically computed critical action 
differences as a function of the pertinent critical energies 
for q — 1. The leading-order result (|26(1 is indicated by the 
solid line. It is accurate even at fairly high energies. 



4 Uniform approximation on the sphere 

In the hydrogen atom in an electric field, the bifurcation 
scenarios are largely shaped by the rotational symmetry 



to? 0.15 

< 




— crit 



Fig. 3. Numerically calculated critical action differences 
(crosses) for q = I — k — lasa function of the critical en- 
ergies. The solid line gives the leading-order result 1261 . 



around the electric field axis: Due to that symmetry, all 
non-axial orbits occur in one-parameter families, whereas 
away from a bifurcation the axial orbits are isolated. A 
normal form describing the bifurcation must also possess 
a rotational symmetry. Therefore, its configuration space 
must be at least two-dimensional. It can be chosen to be 
defined on a Euclidean plane (x, y) and reads [17111811191 



*a(x,») = \r A 



1 2 

—ar 
2 



with a real parameter a and r 2 = x 2 



(27) 

the distance 

from the symmetry centre. Apart from the trivial sta- 
tionary point at the origin, (|27|l has a ring of stationary 
points at r = \/a, which is real if a > and imaginary 
if a < 0. Thus, (|27() describes a family of non-axial real 
orbits present for a > which then contracts onto an 
axial orbit and becomes a family of ghost orbits. It there- 
fore reproduces the elementary bifurcation described in 
section 

In terms of the actions Sax, S non and semiclassical am- 
plitudes vAax, A n on of the axial and non-axial orbits, re- 
spectively, wc found the uniform approximation © corre- 
sponding to the normal form l|27() to be 



9(E) 
where 
and 



^non T 4- — f A -L - — 1 A 
/-< , -\ ~r I Q"/^ax T" r- — -^noii 

(1 + 1) a \ V27T 



i2\/iS' ax Sn 



74 



-c - 



/2vr 



-iS - 



'2tt 



with the Fresnel integrals |2U 
C(x) = j cos(|t 2 ) dt , S(x) 



/ sm — 

Jo V2 



(28) 
(29) 

(30) 



^ dt. 
(31) 
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We saw in section 01 that it is essential at low energies 
to construct a uniform approximation capable of simulta- 
neously describing both bifurcations of a non-axial orbit. 
As the normal form 127(1 describes only one axial orbit 
and a non-axial orbit bifurcating out of it, the uniform 
approximation l|28() is unsuitable for this purpose. 

A normal form describing both the generation and the 
destruction of a non-axial orbit must possess two isolated 
stationary points corresponding to the uphill and down- 
hill orbits and a ring of stationary points related to the 
family of non-axial orbits. As the normal form parameters 
are varied, the ring must branch off the first isolated sta- 
tionary point and later contract onto the second. This is 
impossible if the normal form is to be defined in a plane, 
but it is easily achieved if the normal form is based on 
a sphere. It can then be chosen in such a way that it 
has isolated stationary points at the poles and a ring of 
stationary points moving from one pole to the other. A 
normal form satisfying these requirements is 



(32) 



It is given in terms of the polar angle i? on the sphere. 
Due to the rotational symmetry, it is independent of the 
azimuth angle ip, and it contains two real parameters a 
and b needed to match the two action differences between 
the three closed orbits. 

If, in the vicinity of the poles, <P a ^ is expanded in terms 
of the distance p = sin$ from the polar axis, it reads 



a- 1 



+ _ )+0(p e ) (33) 



around $ = 0, and 

= b(a + l) 2 - ba Qp 4 + ^P 2 ) + O (p 6 ) (34) 

around $ = ir, so that in both cases 1(32(1 reproduces the 
simpler normal form (|27fl up to re-scaling. From the coef- 
ficients of the quadratic terms in 1(33(1 and ((34(1 it is appar- 
ent that a ring of stationary points bifurcates off $ = at 
a = l, and off 1? = 7T at a = — 1. 

These conclusions are confirmed by a discussion of the 
full normal form ((32(1 . The ring of stationary points is 
located at the polar angle i? noll satisfying 



COS $non — 



(35) 



This angle qualitatively corresponds to the starting and 
returning angle of the non-axial orbit in that it describes 
its motion from the downhill orbit at i) — to the uphill 
orbit at ■& = n. Quantitatively, $ non cannot directly be 
identified with the classical angle i? because the energy- 
dependence of the normal form parameters a and b is fixed 
by matching the action differences rather than the angles. 
If —1 < a < 1, the angle i9 non defined by l(3*51) is real. If 
a > 1 , i? non is purely imaginary, corresponding to a ghost 
orbit at energies below E gcn . For a < — 1, # non is of the 
form $aon = 7r + ia, as is characteristic of a ghost orbit at 
energies higher than -Edest- 
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Fig. 4. Bifurcation scenario described by the normal form 132H 
for 6 = 1. Location $ of the stationary points, stationary val- 
ues of the normal form <P at t and stationary values of its sec- 
ond derivative. Solid line: stationary point at l d TLOn , long-dashed 
line: stationary point at ■& = 7T, short-dashed line: stationary 
point at i? = 0. 
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Fig. 5. Orbital parameters for the bifurcations of the (20,21) 
non-axial orbit. Solid line: non-axial orbit (20, 21), long-dashed 
line: 21st repetition of the uphill orbit, short-dashed line: 20th 
repetition of the downhill orbit. The vertical lines indicate the 
bifurcation energies E gcn = —4.155542 and Bdcst = —3.693431, 
respectively. 



The bifurcation scenario described by the normal form 
((32(1 is summarized in figure 0] For 6=1, the stationary 
values of 1(32(1 . the complex polar angles d where station- 
ary points occur, and the second derivative of <P a .b with 
respect to d are shown as functions of a. Due to the ro- 
tational symmetry, the latter is well defined even at the 
poles. 

Figure 0] should be compared to figure [SJ which dis- 
plays the orbital data pertinent to an actual bifurcation 
scenario. The qualitative agreement between figure 0] and 
figurc[S]is obvious, so that ((32(1 can be seen to describe the 
bifurcations correctly. (Note that a is a decreasing func- 
tion of E.) 

As the normal form ((32() has been chosen such as to 
vanish at the non-axial stationary point i9 non , the perti- 
nent uniform approximation reads 



W(E) = I{a,b)e 



iS ao „(E) 



(36) 



where S non (E) denotes the action of the non-axial orbit 
and 



I{a, b) = dd dp P ^)e i ^- bW 



(37) 
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with an amplitude function to be specified below. 

The connection with closed-orbit theory is established 
with the help of a stationary-phase approximation to the 
integral (|37|l . If \a\ > 1, the only real stationary points of 
<P a ,b are located at the poles, and the asymptotic approx- 
imation to i|36|) reads 



6(1 - a) 



6(1 + a) 
If \a\ < 1, the contribution 



P(0)e 
p(ir)e 



i6(l-a) 2 +iS no „ 



i6(l+a) 2 +iS„ 



(38) 



(39) 



has to be added to JSSJ. These asymptotic expressions 
should reproduce the simple closed-orbit theory formulae 



and 



(40) 
(41) 

(42) 

connecting the actions and amplitudes of the downhill, 
uphill and non-axial orbits to the normal form parameters 
a and 6 and to the values of the amplitude function p($). 

Equation l|42[l can be solved for the normal form pa- 
rameters to yield 



A e iS °°° 

respectively. We thus obtain the conditions 

^Sdown = *5down Snon = b(l a) 
AS" up = S up - S'non = 6(1 + a) 2 , 

and 

•^down 



6(1 - a) 



p(0) , A P 



6(1 + a) 



A non = 2lT\l — p(^non) 



down 



and 



clow 

own 



{I -a) 2 (1 + a) 5 



(44) 



(45) 



where the upper or lower sign has to be chosen if the non- 
axial orbit is real or complex, respectively. 

If we set p(-d) equal to a constant in (|43|l . we obtain 
the constraints 



l + i/6 

~. \j "Z «4non — 6(1 dj Adown 

Z V Zn 



1 +i / 6 



(46) 



-6(1 + a) A U r 



the amplitudes must satisfy close to the individual bifur- 
cations. Numerically, l|46fl can be verified to hold in the 
neighbourhood of the bifurcations. 

To obtain a complete uniform approximation, we make 
the ansatz 

p(-&) = pa+pi (cos $ - a) + p 2 (cos 1} - a) 2 (47) 
for the amplitude function. The coefficients are given by 

P0 = P(^non) , 

(1 + a) 2 p(0) - 4qp(t? non ) - (1 - afp^) 

2(1 - a 2 ) ' (48) 

(1 + q)p(O) - 2p(fl non ) + (1 - a)p{n) 
2(1 -a 2 ) 



Pi 



P2 



with the values p(0), p{k) and p(i? n on) determined from 
(I43|l . Due to i|46[) . all coefficients are finite everywhere. 
With this choice, the uniform approximation H36(l reads 



&(E) = p Io +Pih + P2I2 
with the integrals 

f-271 



(49) 



/() 

= 2tT< 



dip / di?sini9e ifc(cos,5 - Q) 
Jo 



C 



C 



-{I -a) +iS 



■(1-a) 



2h M 

-1 + a 

7T 



iS 



2b, 

-1 + a 

7T 



(50) 



in terms of the Fresnel integrals (|3ip. 

f27T 



and 



'0 
7ri 



2;r 



oV / dtfsini?(cosi?-a)e ib(c ° s1? - a) 

p ifc(t+a) 2 _ p i6(l-a) 2 



(51) 



7 2 = / aV / ^sini?(cos^-a) 2 e ib(cos ^ a) 
■/o Jo 

= ((1 - a) e ib ^ 2 + (1 + a) e ib ( 1+a ) 2 ) + ±J . 

(52) 

All integrals are finite as they stand, they do not require 
any regularization. This somewhat atypical behaviour can 
be traced back to the fact that the domain of integration 
in this case is compact. 

As an example, figure compares the uniform approx- 
imation (|49|l for the orbit (20,21) and three different elec- 
tric field strengths to the isolated-orbits approximation 
and to the uniform approximation (|28|l for isolated bifur- 
cations. For the latter, the non-axial orbit and the axial 
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Fig. 6. Uniform approximation for the generation and destruc- 
tion of the (20,21) non-axial orbit for the electric field strengths 
(a) F = 1CT 2 , (b) F = 1(T 6 , (c) F = 1(T 10 . Solid lines: 
Uniform approximation 1491 . long-dashed lines: uniform ap- 
proximation l|28|l for isolated bifurcations, short-dashed lines: 
isolated-orbits approximation. 

orbit for which the action difference is smaller are com- 
bined into the uniform approximation, whereas the second 
axial orbit is treated as isolated. 

For the high electric field strength F = 1CU 2 in fig- 
ure HJa), in the energy range shown none of the uni- 
form approximations reaches the asymptotic region where 
they agree with the isolated-orbits approximation on ei- 
ther side of the bifurcations. Between the two bifurcations, 
the isolated-bifurcations approximation appears just to ar- 
rive at its asymptotic behaviour before the second bifurca- 
tion is felt. The behaviour of the uniform approximation 
on the sphere is completely different there. 

If the field strength is decreased, both uniform approx- 
imations reach their asymptotic regions faster. For F = 
10~ 6 , they approach the isolated-orbits approximation at 



roughly equal pace on both sides of the bifurcations. Be- 
tween the bifurcations, the isolated-bifurcations approxi- 
mation now clearly attains its asymptotic behaviour, where- 
as the uniform approximation on the sphere exhibits qual- 
itatively similar behaviour, but differs significantly from 
the other two in quantitative terms. It requires the even 
lower field strength of F = 1CU 10 used in figure Etc) un- 
til both uniform approximations agree with the isolated- 
orbits result in the region between the bifurcations. 

The uniform approximation on the sphere was con- 
structed because in the energy region between the bifur- 
cations the isolated-bifurcations approximation must fail 
for low energies and high field strengths. However, it turns 
out to reach its asymptotic agreement with the isolated- 
orbits results considerably more slowly than the isolated- 
bifurcations approximation. In cases where the latter re- 
produces the isolated-orbits results between the bifurca- 
tions, of course, no further uniformization is required. 



5 Semiclassical Stark spectra 

More information on the virtue of the uniform approxi- 
mation on the sphere can be obtained from a calculation 
of actual semiclassical spectra. To this end we calculated 
semiclassical photo-absorption spectra for the hydrogen 
atom in an electric field based on closed-orbit theory and 
including uniform approximations. To low-resolution, a 
semiclassical spectrum can be obtained if the contribu- 
tions of orbits up to a maximum period r max are included 
in the closed-orbit sum J5J). 

We will present results for the photoexcitation of the 
Is ground state of the hydrogen atom in an electric field of 
strength F = 1CU 8 a.u.=51.4 V/cm. The exciting photon 
is assumed to be linearly polarized along the field axis. 
In the closed-orbit summation, the uniform approxima- 
tion (|49|l based on the sphere will be used if the action 
differences between a non-axial orbit and both the down- 
hill and uphill orbits it bifurcates out of are smaller than 
2tt. In the other cases, we revert to the simple uniform 
approximation (|28J) or the isolated-orbits approximation, 
as appropriate. 

Figured shows a comparison between the semiclassical 
spectra obtained from the simple uniform approximation 
only (figure [3a)), the closed-orbit summation with the 
uniform approximation based on the sphere (figure Efb)), 
and an appropriately smoothed quantum spectrum (fig- 
ure Efc)). Note that only the oscillatory part of the semi- 
classical spectra is plotted, so that the base line is slightly 
shifted away from zero. This shift is absent in the quantum 
spectrum. 

It has already been found in |21| that, if only the sim- 
ple uniform approximation @ is used, the spectral lines 
close to the centres of the rt-manifolds tend to have con- 
siderable strengths although they are negligibly weak in 
the quantum spectrum. If the uniform approximation on 
the sphere is used, on the one hand the line strengths 
develop a second maximum at the centre of a manifold. 
On the other hand, however, the absolute intensity of the 
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7. Low-resolution photo-absorption spectrum: (a) semi- 
classical spectrum not using and (b) using the uniform approx- 
imation on the sphere, (c) smoothed quantum spectrum. The 
scaled cut-off time is T max = 30. 
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Fig. 8. The manifold n — 23 from the spectra shown in fig- 
ure [7] (a) Semiclassical spectrum not using and (b) using the 
uniform approximation on the sphere, (c) smoothed quantum 
spectrum. 



central lines is considerably reduced, so that the overall 
intensity distribution within an n-manifold is reproduced 
much better than with the isolated-bifurcations approxi- 
mation only. 

These observations can be confirmed if one of the n- 
manifolds is viewed in more detail. Figure |S] compares 
the low-resolution semiclassical spectra of the manifold 
n = 23, which is the lowest n-manifold included in the 
calculation, to the quantum spectrum. In the spectrum 
shown in figure |Bt&), which uses the simple uniform ap- 
proximation (|28|) only, the lines at the edges of the mani- 
fold are the strongest, and the line strengths decrease gen- 



tly toward the centre. The line shapes are slightly asym- 
metric, although the asymmetry is weaker than at higher 
energies ll'll . In the spectrum obtained using the uni- 
form approximation on the sphere (figure|H{b)), the strong 
outer lines of the manifold are virtually identical to those 
in the previous spectrum, but the central lines are dis- 
torted. On the other hand, again, as the central lines are 
weaker in the spectrum of figure IHfb) than in figure [SJa), 
the distribution of intensities within the manifold is de- 
scribed more accurately if the uniform approximation on 
the sphere is used. 
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In addition to the calculation of low-resolution semi- 
classical spectra, we can use the techniques of [2o"ll21] to 
compute individual spectral lines. The results do not de- 
pend strongly on the choice of the uniform approxima- 
tions; the high-resolution spectra obtained from the spec- 
trum in figure [TJb) which use the uniform approximation 
on the sphere (see jHI) are virtually identical to those pre- 
sented in |21) that apply the simple uniform approxima- 
tion (figure [7f a)) only. 

The uniform approximation on the sphere, although 
it satisfies the correct asymptotic conditions and is well 
adapted to the bifurcation scenario qualitatively, in quan- 
titative terms does not yield absolutely precise results. 
The problem is probably connected to the non-trivial spher- 
ical topology of the configuration space of the uniform ap- 
proximation: The construction of uniform approximations 
of the type relies on the assumption that the unknown 
phase function in the diffraction integral can be mapped 
to a certain standard form by a suitable coordinate trans- 
formation. The fundamental theorems of catastrophe the- 
ory 28,29] guarantee that this is possible locally. In a 
Cartesian configuration space and close to a bifurcation, 
all stationary points of the normal form are close to the 
origin, so that coordinate regions away from the origin 
can, in the spirit of the stationary-phase approximation, 
be assumed not to contribute to the diffraction integral. 
A local mapping of the function in the exponent to the 
normal form and a local approximation to the amplitude 
function by means of a Taylor series thus suffice to ob- 
tain good results. By contrast, on the spherical configura- 
tion space a global approximation to the unknown func- 
tions must be sought because the stationary points are 
distributed across the sphere: There is a stationary point 
at each pole. In the case at hand, we tried replacing the 
(cosi?— a) 2 -term in the amplitude function 1)47(1 with sin$, 
which also gives a rough approximation to the actual am- 
plitude function. The alteration has no significant effect 
on the result. The construction of a quantitatively more 
reliable uniform approximation will presumably require 
a deeper global understanding of the bifurcation scenar- 
ios. It remains a worthwhile challenge for future research 
to devise an easy-to-apply technique for the construction 
of uniform approximations that takes full account of the 
classical phase space topology. 

6 Conclusion 

We have discussed the bifurcations of closed orbits in the 
hydrogen atom in an electric field. At low scaled energies 
a uniform semiclassical approximation is required that si- 
multaneously describes both the generation and the de- 
struction of a non-axial orbit. Its construction was carried 
out. Beyond its significance to the Stark system, it is of 
interest for reasons of principle because it uses a sphere 
instead of a Euclidean plane as its configuration space and 
thus is the first uniform approximation introduced in the 
literature that relies on a topologically non-trivial con- 
figuration space. It could be demonstrated that uniform 
approximations with a non-Euclidean configuration space 



can indeed be constructed and allow for the semiclassi- 
cal calculation of significantly improved photo-absorption 
spectra of the hydrogen atom in an electric field. 



References 

1. M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics, Springer- Verlag, New York, 1990. 

2. M. Brack and R. K. Bhaduri, Semiclassical Physics, 
Addison-Wesley Publishing Company, Reading, Mass., 
1997. 

3. M. L. Du and J. B. Delos, Phys. Rev. A 38, 1896 and 1913 
(1988). 

4. E. B. Bogomolny, Sov. Phys. JETP 69, 275 (1989). 

5. A. M. Ozorio de Almeida and J. H. Hannay, J. Phys. A 
20, 5873 (1987). 

6. M. Sieber, J. Phys. A 29, 4715 (1996). 

7. H. Schomerus and M. Sieber, J. Phys. A 30, 4537 (1997). 

8. M. Sieber and H. Schomerus, J. Phys. A 31, 165 (1998). 

9. T. Bartsch, The hydrogen atom in an electric field and 
in crossed electric and magnetic fields: Closed-orbit the- 
ory and semiclassical quantization., Cuvillier, Gottingen, 
Germany, 2002. 

10. T. Bartsch, J. Main, and G. Wunner, Phys. Rev. A, sub- 
mitted. LANL preprint [nTin.CD/02120lTl 

11. J. Main and G. Wunner, Phys. Rev. A 55, 1743 (1997). 

12. H. Schomerus, Europhys. Lett. 38, 423 (1997). 

13. J. Main and G. Wunner, Phys. Rev. E 57, 7325 (1998). 

14. H. Schomerus, J. Phys. A 31, 4167 (1998). 

15. V. Kondratovich and J. B. Delos, Phys. Rev. A 56, R5 

(1997) . 

16. V. Kondratovich and J. B. Delos, Phys. Rev. A 57, 4654 

(1998) . 

17. J. Gao and J. B. Delos, Phys. Rev. A 56, 356 (1997). 

18. J. A. Shaw and F. Robicheaux, Phys. Rev. A 58, 1910 
(1998). 

19. J. A. Shaw and F. Robicheaux, Phys. Rev. A 58, 3561 
(1998). 

20. T. Bartsch, J. Main, and G. Wunner, Phys. Rev. A 66, 
033404 (2002). 

21. T. Bartsch, J. Main, and G. Wunner, J. Phys. B, submit- 
ted. LANL preprint nlin.CD/0212046 

22. J. Gao and J. B. Delos, Phys. Rev. A 49, 869 (1994). 

23. M. Born, Vorlesungen iiber Atommechanik, Springer, 
Berlin, 1925. 

24. L. D. Landau and E. M. Lifshitz, Mechanics, Pergamon 
Press, Oxford, 1960. 

25. M. Abramowitz and I. A. Stegun, Pocketbook of Mathe- 
matical Functions, Verlag Harri Deutsch, Frankfurt/Main, 
1984. 

26. D. Wintgen, J. Phys. B 20, L511 (1987). 

27. M. Brack, Found. Phys. 31, 209 (2001). 

28. D. P. L. Castrigiano and S. A. Hayes, Catastrophe Theory, 
Addison-Wesley Publishing Company, Reading, MA, 1993. 

29. T. Poston and I. Stewart, Catastrophe Theory and its Ap- 
plications, Pitman, Boston, 1978. 



